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I Abstract 

^ I An averaging method for getting uniformly valid asymptotic approximations of the 

' solution of hyperbolic systems of equations is presented. The averaged system of equa- 

tions disintegrates into independent equations for non-resonance systems. We consider 
' the resonance conditions for some classes of solutions. The averaged system can be 

CO , solved numerically in the resonance case. The shallow water problem is considered as 

an example of the resonance system. Results of numerical experiments are presented. 

^ ■ 1 Introduction 

o ■ 

<^ ■ In this paper we consider a system of weakly nonlinear equations with a small positive 

' parameter e: 

O: Ut + A{U)U^=eB{t,x,et,ex,U,U^,U^^„U^^^), (1.1) 



00 



where U = (ui, n2, . . . , Un)"^ is a column vector, A = \\aij\\ is an n x n matrix, and 
bn) is a column vector. We assume that all coefficients are sufficiently 
[ smooth functions. 

^ Many physical problems are described by such systems. We mention only some exam- 

ples: dispersive waves in plasma, problems of weakly nonlinear optics, one-dimensional 
gas dynamic equations, shallow water waves. The construction of uniform asymptotic 
approximations of the solution of (|l.l| ) becomes a nontrivial task if secular terms arise in 
;J] ■ the expansion. We can illustrate this situation by considering the following linear problem 

ut + u^ = eu, 

u{0, x) = sinx. (1-2) 



- 

X 



We find easily the exact solution of problem 

u{t, X, e) = e^* sin(x — t). 
Using a Taylor series expansion of the function e^* we get that 



u{t, X, e) = sin{x — t) + et sin(x — t) -\ — — sm{x — t) + ■ ■ ■ . (1-3) 
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This expansion has secular terms et, e^t^ , . . ., and therefore the formula (|l.3| ) is asymp- 
totical only if et ^ 1. In this case we have that 

u{t, X, e) = sin(x — t) + 0(et). 

For et = 0(1) this expansion is not asymptotical. On the other hand the problem ( |1.2D 
has a classical solution in any domain 

< t + < co/e, 

here cq is a constant. It is easy to prove that 

lim n(t, X, e) ^ sin(x — t). 

Therefore it is a nontrivial task to construct an asymptotical approximation, which is 
uniformly valid in the region t + |x| = O 

The basic idea of all asymptotic methods is to introduce new "slow" variables, e.g., 
r = et, ^ = ex, and to define explicitly the dependence on "fast" variables. For example 



in [21 1 the solution is obtained in the following form: 

Uj = ^j{C, + 0(e), C = r]j = e'^ (x - Xjt + e^^^^^it, x)) . 

Substituting these expressions into system ( |l . l]) , using a Taylor expansion with respect to 
e, and equating coefficients at e^ we get equations for new unknown functions ipj 

Our goal is to approximate the initial nonlinear problem by some fundamental equa- 
tions, for which the analytical solution exists. The Burgers' equation and the Korteweg-de 
Vries equations are examples of such problems. Most asymptotical methods are based on 
physical assumptions and no strict mathematical proofs of the validity of these approx- 
imations are given ||2|, Important applications of asymptotic methods are given in 
1^, |8|, |l9|. A survey of mathematical results on asymptotic expansion methods is pre- 



sented by Kalyakin in [12|, see also Keworkian and Cole [13|. Similar problems arise 
not only for systems of equations, but also for semilinear perturbed wave equations, the 
telegraph equation and weakly nonlinear beam equation 0, 18|. 



Usually application of formal asymptotic method reduces the initial problem to a single 
nonlinear wave equation [ 21 1 . Our method reduces problem ( |l.l[) to the integro-differential 



system of averaged equations. A new problem can be seen as a more difficult problem, 
than the initial formulation. However, we will show that the averaged system disintegrates 
into independent equations in the non-resonance case. Moreover, our method allows to 
describe accurately the resonance interaction of waves, as well. Our internal averaging 
method is close to the multiple-scale method, which is used for similar problems in Q. The 
application of internal averaging method for solving gas dynamic equations is presented 
ini. 

The rest of the paper is organized as follows. In Section ^, we formulate the general 
averaging scheme. We propose the averaging operators and compare our scheme with the 
other similar methods. In Section |3|, we consider the resonance case, when the averaged 
system is connected. Section ^ deals with the application of the proposed algorithm for a 
shallow water model. It is proved that the initial system is ill-posed, hence the regularized 
model is formulated. The averaged system also gives a nontrivial regularization of the 
shallow water model. Finally we give the results of numerical experiments. 



460 



A Krylovas and R Ciegis 



2 Method of averaging 

Let Uq be a constant solution of equation (|1.1[) : 

B{t, x,et, ex, Uo, 0,0,0) = 0. 



(2.1) 



We assume that problem ( |l.l| ) is hyperbolic in the neighborhood of Uq, i.e. there exists 
an n X n matrix R= \\rij\\ det R ^ such that 



A = diag{Xi,X2,...,Xn} = RA{Uo)R-\ 
We are interested in finding a small-amplitude wave solution 

U{t, X, e) = Uq + eUi{t, x, e). 
Linearizing equations ( |L1| ) with respect to Uq yields the system of equations 

Vt + AVz = eF{t, X, et, ex, V, V^, V^^, V^^^) + o{e), 
where V = RUi, and the function F is given by 



(2.2) 



(2.3) 



(2.4) 



F = -RAi [R'^V] R-^V^ + R {Bq [R'^V] + Bi [R'^V^] 

+ B2 [R-^V^cc] + ^3 [R-^V,xx]) , 

d_ 

_duk 



dU 



E 



Bm\Ul] — {bml,bm2, ■ ■ ■ ,b. 

dbj{t, X, et, ex, Uq, 0, 0, 0) 



k=l 

T 

mn j 1 



m = 0,1,2,3 

n 



E 

k=l 



duQk 



Ulk, 



9bj 



k=l 



du, 



Okx 



'Ulkx: 



_ dbj 

Ulkxx, "Sj — / ^ 



k=l " k=l 

If £ = 0, then system ( |2.4D with the initial condition 

V{0,x,e) = Voix) 



duokx 



-Ulkx 



(2.5) 



describes n independent linear waves Vj = voj{x — Xjt). 

If ~ e~^, then the exact solution of initial- value problem (2^), (|2.5D (and therefore 
also of problem ( [LI]) ) is not close to the simple wave. For example, the nonlinear equation 



vt + Vx = evvx 

describes a nonlinear wave, which is given by the implicit relation 

v{t, X, e) = vo{x — t + etv{t, x)). 
Obviously it can not be approximated by a simple wave v = vo{x — t) for et = 0(1). 
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We first rewrite system ( |2.4|) , ( |2.5D in a coordinate form 



dt dx 



efj ( t, X, et, ex, V, 



* ' ' ' ' ' (?x ' dx"^ ' dx'^ J 
Vj{{),x,e) =VQj{ex,x), j = l,2,...,n. 



(2.6) 



Let T = et, = ex be "slow" variables, and yj = x — Xjt, j = 1,2, . . . ,n he "fast" 
characteristic variables. The operator of averaging along the j-th characteristic of the 
non-perturbed (i.e., e = 0) system (2.6) is given by 



Mj[g{t,x,T,^,vi,V2,...,Vn)]= lim. ^ [ g(s,yj + XjS,T,^, 
'"liVj + (^j - Ai)s), . . . ,Vniyj + (Aj - Xn)s)^ds. 



(2.7) 



We define the following averaged system of equations 



dwj dwj 



dr 



. V 9yk dyl dyl J 

n. 



Wj{o,tyj) = voji(,yj), j = i,2 

If the operator B = 0, then averaged system ( ^.8]) takes the form 

dwm 



n n 
^ ^ fjkmMj 
k=l m=l 



Wk 



dyn 



(2i 



(2.9) 



Without a loss of generality we can assume, that 

1 

lim — - / voi((,,x)dx = 0. 
It is easy to prove that the following properties of the averaging operators Mj are valid: 



(2.10) 



Mi 



dvj_ 
^dy- 



Mi 



3 



Mi 



0, 



Vi 



dvj_ 
i +3 



(2.11) 



For n = 2 using ( p.9D -( 2.11 ) we get that the averaged system of equations reduces to 
two independent problems: 



dwj_ 



Wj{Q,i,yj) = voj{^,yj). 



(2.12) 



If the operator i? is a linear function of Uxx or Uxxx, then averaged equations ( p.l2[) 
are described by the Burgers' or the Korteweg-de Vries equations. 

In most cases equations of averaged system (2^) are connected and they describe the 
interaction of waves. Asymptotic analysis is finished at this stage, but some numerical 
analysis is still needed in order to get the solution. 
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3 Approximation accuracy analysis 



Let A4 he a class of functions, for which averages ( |2.7D exist uniformly for all variables. 
Then the solution of averaged system (^]^) exists and we can analyze the approximation 
properties of this solution. 

The case of periodical initial conditions uoj{x) G 6*2^ (M) was considered in [0, For 
sufficiently smooth functions fj{ui, . . . ,Un, uix, • • • , Unx) it was proved that if {vi,V2, ■ ■ ■, 
Vn) is the solution of system (2^) and {wi,W2, . . . ,Wn) is the solution of averaged system 
(|2.8|) , then there exists a constant cq > such that 

limmax sup \uj{t,x,e) — Vj{et,ex,x — Xjt)\ = 0. (3.1) 

I I I — £ 

This result is analogous to the First Bogoliubov's theorem for ordinary differential equa- 
tions [§]. 

The essence of the proposed averaging method is the following: the average Mj[g{T, 
t, X, v)] along the j-th characteristic yj = x — Xjt is the limit of the integral 

1 

hm - / g{T, ^, s, Vj + Xjs, v{t, ^, s, yj + Xjs))ds. 

We integrate function v, which is still unknown, that is why our method can be called 
the internal averaging. For comparison the classical averaging along characteristics [^] is 
defined by: 

1 

hm - / giT,C,s,yj + XjS,v{T,^,t,x))ds 

T-»oo 1 Jq 

and it can be called the external averaging. 

We illustrate the difference between these two techniques, by considering the model 
system with the internal resonance 

ut + Ux = ev sin x, n(0, x, e) = 0, 
vt = 0, u(0, x) = sinx. 

The external averaging gives the problem 

Ut + Ux = Q, [/(0,x) = 

and [7 = does not approximate the exact solution 

u{t, X, e) = ^ (2t + sin 2{x — t) — sin 2x) 

if t ~ e^^. The internal averaging of this system gives the problem 

Ur=lim^[ V{T,y + s)sm{y + s)ds, U{0,y)=0, 
Vr = 0, y(0,x)=sinx, y = X — t. 
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After simple computations we get 



V{t, x) = sinx, 



If 1 

Ur = lim. — sin^ {y + s)ds = -. 

T-*oo 1 Jq I 



Hence U = ^ and the equality u = U + o(l) is satisfied uniformly for t £ [O, O • 

Using Fourier series we can write the resonance conditions for system ( |2.6D . Let 
fj{t, X, r, ^, ui, . . . Vn) be periodic functions with respect to t: 

fj {t + A*, X + Aj^, r, e, vi{t, e, yi + Ai), Vn{T, yn + K)) 
= fj{t,x,T,^,vi{T,^,yi), . . . ,i;„(r,^,y„)). 

Integrating these functions along characteristics we get that 

fj{t,X,T,e,,Vi{T,i,yi), . . . ,Vn{T,^,yn)) 



{lt,l^,h,...,l„)ytO 



2ni 



Then using the substitution 

t = X, x = yj + XjS, yi = yj + {Xj - Xi)s 

we obtain that the condition 

rT 



T 

is satisfied if and only if 



1 P 

hm - / fj{s, yj + XjS, r, ^, . . . yj + {Xj - Xi)s . . .)ds = (3.2) 



4 + f + E^'^^o. 

V;*,r,/fcEZ & \l^\ + \l''\+^\lk\ ^0. (3.3) 



The condition ( ^.2| ) specifies the absence of the resonance in system (p.6|). If it is satisfied, 
then averaged system (2.8) disintegrates into independent equations, which are similar to 

ill)- 

It was proved in [^], that the class of functions A4, for which the internal averaging 
can be applied, includes almost periodic functions. For example, let assume that 

fj = Pj{r,t,,t,x)uiU2- ■ -Un, Pj = 2^ Pj{T,Oe , 
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Then system ( ^.61) is non-resonance if 

Thus our method treats uniformly both the resonance and non-resonance problems. 
The mathematical justification of the method was given only in the case, when the opera- 
tor B in system ( p..l[ ) does not depend on Uxx and Uxxx- But our results for non-resonance 
systems coincide with the results given by some other methods. Therefore we expect that 
asymptotic solution approximates uniformly the exact solution in more general cases too. 
The numerical analysis of such problems will be given in the next section. 



4 Shallow water waves 

In this section we consider the system of shallow water equations 

Zt + {HU)x = e Q {H^Uxx), - \{HU)xxx - HHx{HU)xx - {ZU)x^ , 

Ut + Zx = -eUUx, (4.1) 

where z = eZ is the water surface level, u = eU is the horizontal velocity of the fluid, H is 
the normalized bottom equation. All variables are normalized by some typical horizontal 
(L,.) and vertical (i?*) sizes: 



2 



-iv* -fj* -tj* \ -f^* / 

where xi is the horizontal coordinate, ti the time, zi the water surface equation, g the 
acceleration due to gravity, Hi{x) the bottom equation. We assume, that 

H = l+eh{x). 

Then we can simplify the first equation of the system: 

Zt + Ux = e (^-hjxxx - {hU)x - {ZU)x^ , 

Ut + Zx = -£UUx. (4.2) 

We define new functions and v~ , which are related to U and Z in the following way: 

U = — , Z = + . 
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Then problem (|4.2|) can be reduced to the system 



^f" - = -| ( ^ {vtxx - Vxxx) + {Kx) -V')):, 



')'-{v-)')^-{v+-v-){vt-v-)j. (4.3) 

The asymptotic solution V^{T,y^) = v^{t,x,e) + o{l), = x=ft satisfies the averaged 
system 



dV+ 

dV- 1 /I 



dr 2\3 y 
here we have denoted the averaging operators: 

1 



[Iv-^-y-y- + {{H^)V^).)^_ + 3V-VZ^ , (4.4) 



T 

{hix)V+)_ = hm ^ I h{y- - s)V+{t, y" - 2s)ds, 
(/i(x)y-) = hm i r h{y+ + s)V-{T,y+ + 2s)ds. 

In [l| the analysis of the same problem is based on the assumption that U = Z, i.e. 
only one wave is considered. After simple computations we get the following Korteweg-de 
Vries problem for Z: 

3 1 
Zt Zx + —eZZx + —£Zxxx = 0. 

2 D 



We will prove that in the non-resonance case system (4.4) describes two independent 
waves. 

Let assume, that ( |2.10| ) holds for the initial condition, then we get 

T 



[V^)^=lim^ I V^{T,y)dy = 0. 



System (4.3) is non-resonance if the following equalities 

{h{x)V-)+ = 0, {h{x)V^)-=0 (4.5) 
hold. Then averaged system (^^) reduces to two Korteweg-de Vries equations: 

dr 2 dy 6 dy^ 



466 



A Krylovas and R Ciegis 



Let consider the following initial conditions: 

U{0,x)=0, Z{0,x) r^J^ZkC 



k£Z 



We also assume that 

h{x) ~ ^/ifce*^'=^, 

then the non-resonance condition can be written as 

/Xfc/±2zy,, VA;,/eZ & |A:| + |/|^0. (4.6) 

In the periodic case, when the periods of functions h{x) and Z(0, x) are equal to 2it, we 
have that ^ik = = k. Thus, if h{x) ^ const and Z(0, x) / const, then condition ([4.6|) is 
not satisfied and system (f4.3D has a resonance. 



4.1 Finite difference scheme 

We define the space w/j and time u)t- meshes and assume that the space mesh size h and 
time mesh size r are uniform. We denote by = v{t^,yj) a discrete function defined on 
u)hX uJt-. The following common notation of difference derivatives is used in our paper 

V'^'^'^ — Vj — Vj-l 

Vr = - , Vy = - , 

_ Vj+1 - Vj _ Vj+1 - Vj-l 
On — : , (Jo — . 

^ h ' y 2h 

The finite difference approximation of system (^^) is defined as follows (see also 

1 /y'l+i + 3 / (y"+i)^ + + (y" 



F+ {W^+^,W,j + 1) - F+ (PF"+\Ty",i - 1) 

_ 1 + 3 / [W+^f + W+^W + (VF")^\ 

" 6 I ^ + 4 1^ 3 J . 

F_ (y'^+i, + 1) - F_ (y"+^y^ j - 1) 
+ 4/i ' ^ ' 

where the integrals are approximated as follows: 

1 ^ T/^+l -I- T/" 



1=1 



The approximation error of this finite difference scheme is estimated as O (r^ + /i'^). Nu- 
merical methods for solving the Korteweg-de Vries equation are investigated in |^, 17|. 
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4.2 Linear dispersion problem 

In this section we consider a linear problem 

£ 

Zt + {HU)x = — - Uxxx, 

Ut + Zx = 0. (4.8) 

First we will prove that system (^]^) defines an ill-posed problem. Let consider the case 
H = 1. After simple computations we get the equation for U : 

Uu -Uxx = - Uxx... (4.9) 



Considering the k-th. Fourier mode we get that the solution of (4.9) is unstable for k'^e > 3. 
In order to define a stable solution we use the following regularized problem 



Zt ~\~ {HU^x — 2 Uxxx 20 ^XXXXXl 

Ut + Zx = 0. (4.10) 



We note that the averaged system (4.4) also gives a nontrivial regularization of this 
ill-posed problem. 

The accuracy of asymptotic solution is illustrated by solving problem ( |4.10| ) with fol- 
lowing initial conditions 

U{x,0)=0, Z(x, 0) = cos X + sin 2a;, h{x) = 5sm2x. (4-11) 

Fig. |l] shows the solution of system ( [I.IOD and the asymptotic solution at t = 1/e, for 
e = 0.1 and e = 0.01. For e = 0.001 the difference between the asymptotic solution 
and the exact solution is too small to be illustrated on the figure. We also note, that 
the averaged system must be solved numerically only once and then a solution can be 
computed using simple interpolation procedure for any parameters e. 
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4.3 Nonlinear nondispersive problem 

In this section we consider a nonlinear problem 



Ut + Z, = -eUU, 



(4.12) 



with the same initial conditions ( 4.11 ). As it follows from results given in previous sections, 
we have the resonance case. Fig. ^ shows the solution of system ( [4.12| ) and the asymptotic 
solution at t = 1/e. 
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Figure 2. Asymptotical solution for nonlinear problem (4.12). 



4.4 Shallow water waves 

Here we solved the full nonlinear system ( |4.1| ). The regularization of sectio n |4.2| is used 
in order to define the stable solution. Fig. ^ shows the solution of system ( [4.lD and the 
asymptotic solution at t = 1/e. 




Figure 3. Asymptotical solution for the shallow water problem. 
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